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Abstract 

Plane- wave reflection and refraction at an interface with a double wire medium is considered. The 
problem of additional boundary conditions (ABC) in application to wire media is discussed and an ABC- 
free approach, known in the solid state physics, is used. Expressions for the fields and Poynting vectors 
of the refracted waves are derived. Directions and values of the power density flow of the refracted waves 
are found and the conservation of the power flow through the interface is checked. The difference between 
the results, given by the conventional model of wire media and the model, properly taking into account 
spatial dispersion, is discussed. 

PACS numbers: 41.20.Jb, 42.70.Qs, 77.22. Ch, 77.84.Lf 




FIG. 1: Geometry of the unbounded DWM. 



I. INTRODUCTION 

Wire medium (WM) is an artificial medium formed by a lattice of ideally conducting thin wires. 
Recently, we have observed a growing interest to such artificial media of possible new physical 
phenomena and potential applications. This medium at low frequencies is usually described as 
a uniaxial crystal, whose permittivity tensor components are expressed by the plasma model. It 
has been shown |l| , that if the wavevector in a WM has a nonzero component along the wires, the 
plasma model should be corrected introducing spatial dispersion (SD). Similarly, spatial dispersion 
is inherent to the double WM (DWM) formed by two mutually orthogonal lattices of thin ideally 
conducting straight wires, see Fig. [T] 

For consideration of waves in double wire media let us take a case where the wires are per- 
pendicular to each other in the y and z directions. The waves in unbounded space filled with a 
DWM medium were studied in |2fl numerically, in (^J using a semi-analytical approximation of the 
local field, and in ^| both numerically and using the effective medium (EM) approach. In the last 
paper a very good agreement between the results given by the EM and full-wave theories for all 
types of waves in DWM (if the wires are thin) has been demonstrated. 

In this paper we consider the plane-wave reflection and refraction at an interface of DWM using 
the effective medium approach. We assume that the two orthogonal wire arrays are identical, the 
period of the lattice is equal to L in the x, y, and z directions, and the radius of the wires is equal 
to vq. In this case the wire lattice is square in the plane of the wires, i.e., the (yz) plane. We 
assume also that the interface of DWM lies in the (xy) plane and the incident wave vector lies in 
the (yz) plane (see Fig. EJ). 

In the next sections we will give basic expressions, obtained in the framework of the EM ap- 
proach, and formulate the wave refraction problem, demonstrating necessity of additional bound- 
ary conditions (ABC's) both for interface problems with single and double wire media. Then 
we will discuss some approaches, applying in solid state physics in order to overcome the ABC 
problem and look what may be useful for us in application to WM. Using the ABC-free approach 
we will find the reflection coefficient, amplitudes and Poynting vectors for refracted waves. We will 
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FIG. 2: Geometry of the wave reflection problem. 

compare directions of the group velocity and the energy density flow found from the expression 
for the Poynting vector containing additional terms inherent for media with spatial dispersion. 
Conservation of the normal component of the power flow vector at passing through the interface 
is checked. 



II. FIELD EQUATIONS AND EIGENWAVES IN DOUBLE UNBOUNDED WIRE 
MEDIUM 

Assuming space-time dependence of fields as e l ^ Jt ~ kvV ~ kzZ ' , there are non-zero wave vector 
components parallel to wires. Anisotropy appears in this electromagnetic crystal with square 
lattice and DWM behaves as a biaxial crystal with the relative permittivity dyadic 

e = e x u x u x + e y u y u y + e z u z u z , (1) 

with 

e x = e h , e y = e h (l- p-f^|J > e z = e h (l- p"f^|J , (2) 

where k = -y/ch, c is the speed of light, and is the relative permittivity of the host medium. 
Note, that model (EJ) works both for real and imaginary k y (for propagating and evanescent waves, 
respectively), see 

In the wire medium the Maxwell equations 

V x E = -jujfi H (3) 

V x H = jooe f ■ E (4) 

split into two separate subsystems describing wave propagation of fields with TE and TM polar- 
izations: 

-j(k y Uy + k z u z ) x (E x u x + E y \iy + E Z U Z ) = -jufi (H x u x + HyUy + H z u z ) (5) 
-j(kyU v + k z u z ) x (H x u x + Hy\iy + H ' Z U Z ) = j w e Je x E x u x + e v E v u y + e z E z u z ). (6) 



For ordinary (TE) waves this leads to the wave equation 



[k 2 -k\- kl]E x = 0. (7) 

The same equations hold for H y and H z . There are no effects due to wires, and ordinary waves 
propagate as in any isotropic dielectric medium. 

Whereas for extraordinary (TM) waves the wires affect the propagation, and we obtain the 
wave equation 

[k\ - k 2 z - ky—]H x = 0. (8) 

The same equation can be written for E y and E z . 

In order to solve the wave reflection problem we need to evaluate the eigenwaves which are 
outgoing from the interface of DWM. It means that we have to find k z under fixed k and k y . Let 
ky — — k sin 9, where 9 is the incidence angle. 

It follows from (JHJ) that the dispersion equation has the form 



k dy{ky) k. 



2 €y(k y ) 

y e z {k z )\ 



0, (9) 



T(k z ,uj) = k 2 z 
and its solution is: 

2k A - 2k 2 k 2 p - 3k 2 k 2 y + 2k 2 p k 2 y + k 4 y ± kyyjik* - Wyjlkl + k 2 y ) 2 - k 2 (4k 2 + k 2 )) 
k " z± = 2(k 2 - kj) ' (10) 

Two waves propagating or attenuating in both directions follow from the effective medium theory 
(jiup. The conventional isotropic plasma model leads to only two waves for a certain direction, 



namely, k z = ±^jk 2 e — k 2 , e = e^(l — k 2 /k 2 ), where k = uj/c. 

For the following consideration of the refraction problem we will need to know properties of 
waves propagating in the z direction. Here we briefly revise these properties (see our previous 
paper for more details). The real and imaginary parts of k z versus the normalized frequency 
k/k p are presented in Fig. 01 for 9 = tt/4. Let us assume that = 1. In the simple case of the 
conventional model (dotted curve) k z is imaginary if k < K 2 = k p / cos 9, and it is real if k > K 2 . 
The solution of the conventional model is completely wrong between K\ and K 2 because it gives 
an imaginary value of k z instead of a real one which is obtained from the correct model. The 
correct, more complicated solutions, follow from Eq. JTUJ). Analyzing Eq. (|1U|). one can see that 
there exist three frequency regions, corresponding to different kinds of solution. 

The first one is the low frequency band k < K\, where 



K 1 = k p ^-J — 11 

sm9 V cost^ 

is the stop band edge: Beyond this wave number the waves are propagating. There the propagation 
constant k z is complex despite the fact that we have assumed the medium to be lossless (see Fig.|HJ). 
Actually, there are two complex conjugate solutions for each ~Re(k z ) > 0. 




FIG. 3: Real and imaginary parts of k z , calculated using the electro dynamical model (solid curves) and 
the EM theory (dashed curves). The dotted curve shows k z given by the conventional plasma model. 

The second frequency area is K\ < k < K^. At point K2 one of the solutions is zero and 
within the range K\ < k < K 2 we have a forward wave and a backward wave with respect to 
the interface. It means that one of the waves has a positive projection of the wave vector on the 
interface inner normal, while the other wave has a negative projection. 

Finally, for k > K 2 both of the waves are propagating forward waves. Electrodynamical 
calculations 0] (using the three dimensional Green's function) confirm the results of the effective 
medium theory with a high accuracy in a wide spectral range including the regions of evanescent 
and propagating waves (see solid and dashed curves in Fig. EJ). Note that point K2 corresponds 
to the edge of the passband in the framework of the old model. Thus, the model taking into 
account spatial dispersion leads to a considerably more complicated structure of eigenwaves than 
the conventional model of an isotropic plasma, and it is in very good agreement with the results 
of the full-wave analysis. 



III. WAVE REFLECTION FROM A WIRE MEDIUM INTERFACE AND THE PROB- 
LEM OF ADDITIONAL BOUNDARY CONDITIONS 

As it was shown above, there exist two extraordinary waves with the wave vector and the electric 
field in the (yz) plane. Assuming the y-component of the electric field of the incident wave to be 
equal to unity, and applying the continuity conditions of the tangential field components results 
in the reflection problem formulated as follows: 

1 + R E = E + + E_ 

(1 - R E )/Z = E + /Z + + 

where Re is the unknown reflection coefficient for electric field, E + , E_ are the unknown ampli- 
tudes of refracted waves in the wire medium, Zq is the wave impedance (TM) in free space and 



Z± are the wave impedances of the refracted waves. Thus the problem becomes similar to one 
appearing in crystallooptics, where excitons arise and spatial dispersion cannot be neglected Q. 
The main difficulty here is the necessity to invoke additional boundary conditions (ABC) in order 
to match solutions at the interface of media. It was pointed out first by S. Pekar |8j (1956), that 
the well-known Maxwell's boundary conditions (fT2*|) are not sufficient to connect the amplitudes 
of the incident and transmitted waves in adjoining media, if more than one independent wave can 
propagate in any medium. 

Probably, the first ABC were proposed by S. Pekar p], and his ABC stay among the more 
often used in the theory of media with SD. The simplest phenomenological Pekar's ABC are the 
following: 

P = 0, 

S = o (13) 

at the boundary, where P is the polarization vector. Let us consider if we can apply to the WM the 
first condition of (|13|) in the region of evanescent waves. For the ^/-component of the polarization 
vector one can write 

P y = XyiEi + X y2E 2 = 0, (14) 

where the susceptibilities Xyi an d Xy2 relate to different waves, but Xyi — Xyi- Hence we come to 
E\ + E2 — 0, which means that the reflection coefficient for the electric field Re — — 1. In other 
words, a semi-infinite WM behaves as a perfect metal plane. However, we know that it is not so 
even if only evanescent waves are excited in the WM. The reflected wave changes its phase due to 
penetration into WM, and this phase is not equal to 71. The same result is obtained considering 
condition P z = 0. The second Pekar's condition (|TH|) (for derivatives) leads to a similar result. 
Thus, the conventional Pekar's ABC are useless for solving our problem. 

After pioneering Pekar's publication different ABC were proposed for problems of crystalloop- 
tics, as well as semiconductor and plasma electrodynamics. All of these works relate to specific 
media and take into account the properties of a sub-surface layer at the media interface (see BO 
and the bibliography presented there). Besides, phenomenological assumptions and experimental 
data are used in these theories. Since WM strongly differs from a solid-state crystal, only general 
approaches to ABC which do not concerned with particular media, may be interesting for us. 



Many authors (see, for example, 10), [llj) derived ABC from a given model of medium with 



an explicit specification of its surface. After simplifications these models give nothing more than 
Pekar's ABC. Among different theories of SD media so-called 'ABC-Free" theories attract our 
attention. One of the ways for ABC derivation is a concept of "exciton dead layer", proposed by 
Hopfield and Thomas. Actually, there is a layer in which exciton wave function has evanescent 
form. Kikuo Cho 12 1 has pointed out that there is a certain limit of transition layer thickness, 



below which the ABC theory is unnecessary. However, Cho theory is not yet free from some 
parameters determined by specific surfaces. B. Chen and D.F. Nelson declare that their work 
ft solves the macroscopic ABC prob.em completely. Desp.te that the authors of Q used a 
complicated quantum mechanical model, their derivation found that the fully macroscopic solution 



is equivalent to the use of Pekar's ABC P(z = 0) = 0. Thus, such an approach also is not suitable 
for us. A. Vinogradov and co-authors, considering in effects of SD in composite metamaterials, 
have obtained ABCs using an assumption of existence of additional waves in free space which are 
the same as in the metamaterial but have evanescent nature. This assumption leads to additional 
equations at the interface. 

Another way to solve the problem was proposed recently by K. Henneberger It is based on 
the assumption of an abrupt transition from medium to vacuum. It is assumed that the incident 
wave excites a source s(z,u), located within a sub-surface layer < z < 2a, and its thickness is 
assumed to be negligible. This approach is appropriate for our problem of wire media interface. 
Indeed, it is known for problems of diffraction by single semi-infinite wire grids that the induced 
currents deviate from the regular amplitudes far from the interface only in a very narrow region 
whose width is of the order of the grating period [16J. This conclusion holds for grids of wires 
both parallel and perpendicular to the edge. For this reason we assume that for the wire medium 
interface the transition layer has the thickness of only a few periods of the lattice. This thickness is 
much smaller than the wavelength and negligible as compared with the length of the wires (wires 
are infinite in our model). 

Applying this approach to our interface problem of free space and the wire medium, the wave 
equation for H x in unbounded medium (Eq. (jHj) written for H x ) should be replaced by an inho- 
mogeneous equation 

d 2 H T 



dz 2 



x + 



H x = s(z,u), (15) 



'e z (k z )_ 

where H x is the refracted field. It means that any propagating wave has to be created by a source. 
The proper source of the penetrating wave in the wire medium is the incident wave and the 
polarization induced by it in the medium. Such an externally controlled source can be identified 
with some polarization additionally induced to the one already described by e. Therefore it is 
located only on the surface and in the transition region, where the induced polarization deviates 
from that in the bulk medium. After Fourier transform of Eq. (fTB^) one obtains 

H x {z,uj)= — — , 16 

i_ 0O 27r T(q,w) 

where s(q,u) is the Fourier transform of s(z,u) and T(q,u) is determined by Eq. ©. Assuming 
an abrupt transition from the medium to vacuum, we can present the source as a delta function 
s(z,u) = Sq(u)S(z), then s(q,u) = s (u>). If T(q,ou) is an analytical function, the integration in 
Eq. (dnj) can be performed using the residue method. Residues can be found by presenting 1/T 
in the form 

" + + (17) 



where the coefficients are 



£-2 _ U2 ' U2 _ p 



h 2 — h 2 — h 2 

P _ h 2 ' 

^z- fXj z+ 



The residues (the relative amplitudes of the transmitted field components) read 



h 2 - h 2 - h 2 



2k z + 2{k z _ k z _^)k z 



(19) 



(20) 



P- 



h 2 - h 2 - h 2 



2k z - 2(k z _ k z +)k z — 
Finally, the field component H x in the wire medium is 



H x = s [R+e- jk *+ z + R_e~ jk *- Z ] = s Q 



P + e ~jkz+z _|_ P- Q -jk z -z 



2k 



z+ 



2k 7 



The expression for E y and E z can be obtained from the Maxwell equations as 

k z H x = —ko6 y Ey, —k y H x = —koe z E z , 
which gives us the electric field components in the wire medium 

E y = --^- [P + e~ jk *+ z + p.e~^~ z ] 

ko€y 



(21) 



(22) 



(23) 



(24) 



_ kySp 

uje n 



P, 



P- 



_ e ~i k z+z _|_ f c ~.jk z -z 

2k z _e z _ 



(25) 

_2k z+ e z+ 

Now we have expressions for all field components induced in the wire medium for the TM polar 
ization. 

In free space there exist incident and reflected TM waves. Magnetic field components are 



HI = H e- jk « y e- j P° z , H r x = H r e~^ v ^ z = R H H e- jk y y e j ^ z , 
where p = ^Jk^ — k 2 , and the electric field components are 

Ei = -^H e~ jk y y e- j(3az , E r v = R H ^H e- jk y y e^ oZ 
y h y k 



(26) 



(27) 



El = ^H e- jk v y e- jf3oZ , E r z = R H ^-H e- jk y y e jf3 ° z . 
ko k 

At the interface z = the continuity of the tangential field components leads to relations 



Hq + RrH = y 



P ±+ P z 
k z + k z - 



(28) 



(29) 



(30) 



from which the reflection coefficient for magnetic field is obtained: 

§+ I JL_\ _ i_ (§+ , f± 

k z+ k z — J e<y \ Po Po 

Finally, the explicit expression for the coefficient s (the transmission source) is obtained: 

S ° = +"*=-) H ° = T^TMTTT^T^- (32) 

2^k z+ ^ k z _> T— + T— )+— flT + ^ 



k z+ k z _ J e y y/3 (3 

In the region of complex waves k/k p < K\ we have to chose the branches of the square roots for 
k z ± having positive imaginary parts. In the region K\ < k/k p < K 2 it is necessary to take for the 
backward wave ( " — " wave) the root branch with the opposite sign. In the above derivation we 
evaluated the reflection coefficient for the magnetic field. The reflection coefficient for the electric 
field is Re = —Rh cos 6. 

A. Reflection from a single wire medium interface 

The plane wave reflection coefficient from an interface of a single wire medium where wires 
are along the z axis is easily obtained as a special case of the previously considered double wire 
medium reflection problem. In a single wire medium the permittivity components are e x = e y = 
and e z = e h ^1 — k z\z j ■ Evaluating the dispersion equation we have as solutions k z+ = k, which 
is the propagation factor for the TEM mode and k z _ = ^Jk 2 — k 2 — k 2 for the TM mode. Thus, 
in the single wire medium the two extraordinary eigenwaves are TEM and TM polarized. 

We can use exactly the same expressions for the reflection coefficient as in the case of the double 
wire medium simply substituting k z+ = k and k z _ = ^Jk 2 — k 2 — k 2 . This leads to expressions 
for the coefficients 

^ = ¥TV 2 (33) 

^wk (34) 



and residues 



R+ ~ 2A^+ W) (35) 
k 2 

R. = — y . (36) 

9 /P _ P — P U2 _i_ 1,2^ v ' 



IV. DISCUSSION 



Fig.HJshows the phase of the reflected wave <fi = arctan{Im(i2#)/Re(.R#)} versus the incidence 
angle 9, calculated for different k/k p , corresponding to the region of complex waves (see Fig. EJ). 



Here Re(-R#) and Im(i?//) are the real and imaginary parts of the reflection coefficient Rh- The 
calculations have been performed using the model (J2J), taking into account spatial dispersion, and 
the conventional one, 

e 



(37) 



As it is expected, both models give the modulus of the reflection coefficient equal to unity due to 
the absence of losses and propagating waves in wire medium. It is remarkable, that the spatial 
dispersion results in a weak dependence of the phase of the reflected wave on the incidence angle. 
The parameters of the wire medium are taken as above in the eigenvalues calculations. 

The real and imaginary parts of the reflection coefficient in a wide spectral range cover areas 
with complex waves, FW and BW waves and FW waves only, are presented in Fig. El The incidence 
angle is taken to be 7r/4. For understanding of these characteristics it is useful to compare them 
with the eigenvalue dispersion (Fig. EJ). Distinctive points K±, where propagating waves appear 
in the framework of new model, and K 2 , where both of the waves become forward ones. Results, 
given by the old model, are also shown here. The most important feature is that the area of 
propagating waves shifts by Ak/k p = K 2 — K\ in comparison with that given by the old model. 




FIG. 4: Normalized phase of the reflected wave 4>/ti versus the incidence angle 6 /it, calculated at different 
k/k p : curves 1 correspond to k/k p =0.3 and curves 2 correspond to k/k p =0.7. Solid and dashed curves 
show the reflection phase given by the new and conventional models, respectively. 



V. GROUP VELOCITY, POYNTING VECTOR AND REFRACTED WAVES IN DWM 



In this section we will discuss the group velocity and Poynting vectors of waves excited in a 
double wire medium and check the power conservation at the interface. It is well known that the 




FIG. 5: Real and imaginary parts of the reflection coefficient versus the normalized wavevector k/k p . 
The solid line shows Re(i?£;), the dashed line shows Im(i?^) (the new model), the dash-dot and dot lines 
correspond to Re(i?£) and Im(i?£), respectively, obtained in the framework of the old model. 



group velocity is defined as 



v gr = grad k u;. 



(3* 



The Poynting vector that determines the energy density flow for media with spatial dispersion has 
the form 171 ]: 

1„ ™, ujde ik 



S = -Re{E x H*} - 



-E*E, 



(39) 



2 — i > 4 4 ' 

where the permittivity dyadic components are expressed by formulae (J2J) for DWM. Their partial 
derivatives read 



de x 
dk a 



0. 



d€y 

dk,, 



2k 2 k y 



de, 



2k 2 k z 



(k 2 -k 2 ) 2 ' dk z {k 2 -kf) 2 



(40) 



Let us derive expressions for the Poynting vector in free space and in the wire medium. In free 
space the field expressions are 



(41) 



k,H n 



k k 1 

The fields of the waves, marked by + and — , in the wire medium are 

so Pi 



H 



2± 



2 k 



!± e- jkz±z e~ jkyy u, 



E 



2± 



spy 

2k e v 



jk z ±z f ,-jk y y u + K S 0V P± r -jk g ±z r -jk v y 



2k k z ±e z ± 



(42) 

(43) 
(44) 



Now we can write the Poynting vector in free space and in the wire medium (at the interface): 



S 1 (0) = -Re{E 1 xH*} 



\H \ 2 V 



[A)(l - \Rh\ 2 )u z + ky(l + 2\R H \ cos0 + \R H \ 2 )u y ] (45) 



2~~~ L ' 1J 2k 
with the notation Rh = \Rh\&^- In the wire medium, the cross product term is 



S° ± (0) = -Re{E 2± x H; ± } 



— Re < u z + k v - — : 

2k 4 {e y k* z± k z ±k* z± e z ± 



u, 



(46) 



and the spatial dispersive term is 

\s \ 2 r) 1 



2± 



(0) 



2k 4 



u z + 



kpky 



(U2 _ h.2 \2U , U* e , p* z (U2 _ h^c 2 
H z ±) K z ± K z± e z±£ z ± l ft K y) e y 



;f3±(3±U y 



(47) 



The total Poynting vector in wire medium is 

S 2± (0) = S° ± (0) + S^ ± (0). 



(48 



As an example, we consider excitation of modes, located above the plasma frequency Q. For 
the taken parameters of DWM, the wire radius ro=0.01 cm, and the period L=l cm, the plasma 
wavenumber is k p ~ 1.38 cm -1 . Dispersion diagram in form of isofrequencies and directions of the 
Poynting vector (the energy velocity) for these modes with respect to the normal to the interface 
(z-axis) are shown in Fig. |H1 The calculations were performed at k/k p = 1.35. The value of 
the tangential component k y is determined by the incidence angle 9, namely, k y = ksinO. The 
direction of the group velocity, found by numerical differentiation of the dispersion characteristics, 
exactly coincides with the direction of S2±, calculated using formulae (}4l)|) - (j4*%j) . Disregarding the 
term that takes into account spatial dispersion leads to a strongly incorrect result, illustrated by 
dashed curves. 




o Km 0.2 



FIG. 6: 1. Isofrequencies, dotted curves. 2. Solid curves show the angle ip between the energy velocity 
and the z axis versus k y . 3. Dashed curves show the angle between S2±(0) and the z axis versus k y . 



Next, let us consider the power conservation at the interface of free space and the wire medium. 
It is important to check this properly because in the wire medium region there exist two waves, 
and the wire medium is spatially dispersive. In the frequency range considered here the parameter 



values are assumed to be real. The normal components of the Poynting vector on both sides of 
the interface are (using the continuity condition of the tangential field components) 



and 



5-2,(0) 



lUJ6 r 



S zl (0) 



k y _i_ k-, _ 
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8u;e e, 
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k z + k z . 



kpky(3+ 
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Subtracting these power density components leads to expression 



■Mo) - s z2 {o) 



8ue n 
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h 2 h 2 

It p K y 



(k z _ k z 



(49) 



(50) 



(51) 



Using the expressions for f}±, k 2 z ± and e y , we find that the term inside the square brackets vanishes. 
The normal component of the Poynting vector is continuous across the interface, which means that 
the power conservation law is satisfied. 

The normal to the interface components of the Poynting vectors of the refracted waves as well 
as the reflection coefficient (which is purely real beyond Ki) are shown in Fig. [7| The values of 
S2Z+ and S2z- are normalized to the power density of the incident plane wave \Si\ — ^i]Hq cos 8. 
These results agree with those shown in Fig. El namely, S% z - becomes zero at point k/k p = K 2 . 
The point K 2 is a transition point of the " — " wave because it is a backward wave with respect 
to the interface if k/k p < K2, and a forward wave if k/k p > K2. 
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FIG. 7: Reflection coefficient Rh (dotted curve), £22+ (solid curve) and 5*22- (dashed curve) versus the 
normalized wave vector. 



Angular dependence of the normal components of the Poynting vector is presented in Fig. |HJ 
Both refracted modes are TM modes in this case. At low frequencies (k/k p < 1) a plane wave 
cannot excite propagating TM modes at any incidence angle 8, so we present here the results for 
k/kp = 1.5. There exists a critical angle 8 C ~ 0.37T, such that S z \ = S z2 — if 9 > 8 C . As in the 
previous case, Rh = +1 if 8 = 8 C , and Rh — > 1 if 9 — > tt/2, so DWM is an electric wall for grazing 
incidence. 



FIG. 8: The real and imaginary parts of Rh (dashed and dotted curves, respectively), S z \ and S Z 2 (solid 
curves) versus the incidence angle, calculated at k/k p = 1.5. 



A. Poynting vector in single wire media 

The Poynting vector in a single wire medium can be simply deduced from the expressions for 
the double wire media. In this case the partial derivatives of the permittivity components read 



de x 
dk x 



0. 



d€y 

dk,„ 



0. 



dk z 



2k 2 k z 



(k 2 -k 2 z y 



(52) 



The additional term arises now only in the z component of the Poynting vector, and it is equal to 
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Thus, the total Poynting vector reads for real parameter values 
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This gives us after substituting k z ± and f3~ 
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(55) 



(56) 



It is known that for each eigenwave the group velocity in single wire media is parallel to the 
phase velocity. From the above expressions one can see that the Poynting vector is also in the 
same direction as the group velocity (and the phase velocity). The continuity condition of the 
power flow across the interface (continuity of the normal component of the Poynting vector) can 
be easily checked similarly as for double wire media. 

Fig. El illustrates the angular dependence of the reflection coefficient Rh and the normal to the 
interface component of the Poynting vector S z for the TEM mode (Stem)- Parameters of the WM 
are taken the same as in the case of DWM. The amplitude of the magnetic field of the incident 



FIG. 9: The real and imaginary parts of Rh (dashed and dotted curves, respectively), Stem (solid 
curves) versus the incidence angle, calculated at k/k p = 0.5. 
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FIG. 10: The real and imaginary parts of Rh (dashed and dotted curves, respectively), Stem and S z tm 
(solid curves) versus the incidence angle, calculated at k/k p = 1.5. 

wave is H = 1, and the calculated Poynting vector is normalized to the incident wave power flow 
density. Wavenumber k/k p corresponds to the frequency below the plasma resonance and the TM 
mode cannot be excited at any 9. For grazing incidence Rh — > — 1. It is obvious that in this 
case the surface must behave as a wall, because no traveling waves in the medium can be excited. 
Actually, it behaves as a magnetic wall, because in this geometry the wires are orthogonal to the 
interface, and the normal component of the electric field is zero. 

More complicated is the case when the wavenumber is taken beyond the plasma resonance, i.e. 
k/kp = 1.5, see Fig. El However, it follows from Eq. (j2J), that for oblique wave incidence, the 
"effective" plasma wavenumber becomes k' p = k v j cos 8 because k y = ksin.6. For grazing incidence 
Rh — > — 1, similarly to the previous case. 

For small incidence angles both TEM and TM modes are excited. There exists an incidence 
angle 8', which is the angle of total reflection. It takes place when k = k' p , so cos #'=2/3. Neither 
TEM nor TM modes are excited, and Rh —>■ 1, which corresponds to the case of an electric wall. 
For larger 9 the TM mode disappears. 



VI. CONCLUSIONS 



Considering the problem of a plane wave refraction at the interface of a double wire medium, 
exhibiting strong spatial dispersion, we have shown, that Pekar's additional boundary conditions 
are not applicable for its solution. We have analyzed the literature discussion on the ABC problem 
and have come to the point of view, that Hennenberger's approach 4.5] for SD media can be applied 
for some kinds of metamaterials in the microwave range including single and double wire media. 

Despite of criticism of the used ABC-free method, see 3, n| and also replies 2(| , we suppose 



that application of this approach has allowed us to overcome the problem of additional waves 
appearing in media with spatial dispersion and to obtain reasonable results for reflection coeffi- 
cient and for the power densities and group velocities of the refracted waves. Fulfilment of the 
conservation of the power, passing through the interface, can be considered as an evidence of the 
correctness of this approach. 

Possible applications of DWM include antenna structures, low-frequency filters, frequency se- 
lective radomes, double negative metamaterials. 
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